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Abstract 

A unified integrable system, generating a new series of interacting matter-radiation models with 
I interatomic coupling and different atomic frequencies, is constructed and exactly solved through 

algebraic Bethe ansatz. Novel features in Rabi oscillation and vacuum Rabi splitting are shown 
, on the example of an integrable two-atom Buck-Sukumar model with resolution of some important 

' controversies in the Bethe ansatz solution including its possible degeneracy for such models. 

O ' PACS numbers: 02.30 Ik, 42.50 Pq, 03.65 Fd, 32.80 -t 

. The basic physics underlying a variety of important phenomena in interacting matter-radiation 

, (MR) systems, like those in quantum optics induced by resonance interaction between atom and a 
! quantized laser field, in cavity QED rempe879p,raizen89, in trapped ion interacting with its center 
of mass motion irradiated by a laser beam trap,vpgel95 etc., seems to be nicely captured by simple 
models like Jaynes-Cummings (JC) jp. Buck Sukumar (BS) bs and some of their extensions jql. Many 



C ' theoretical predictions based on these models, like vacuum Rabi splitting (VRS) raizen§9,VRS, Rabi 
^ ' oscillation and its quantum collapse and revival rempg8790 etc. have been verified in maser and laser 
O^' experiments. However, for describing physical situations more accurately one has to look for further 
generalizations of the basic models, like q-deformed BS and JC model qbs^qjc, trapped ion (TI) with 
^ , nonlinear coupling vogel9§,ntrap, multi-atom models raizen89,ntrap,natom etc. Nevertheless, while 
^ I the exact solutions for the JC and the BS models together with their simple multi-atom extensions 
are known jcexact,jcbethe,bsbethe, the same is no longer true for most of the above generalizations. 
Moreover, while in known multi-atomic MR models the atoms interact only via the oscillator mode 
jcbethe^bsbethe with coinciding atomic frequencies (AF), integrable models with explicit inter-atomic 
couplings have not been proposed. Likewise, though q-deformation, which physically signifies intro- 
duction of anisotropy together with specific nonlinearity into the system, was considered for a few 
MR models qbs^qjc, their multi-atom and integrable variants are not known. Therefore, it is indeed a 
challenge to find a scheme for generating integrable MR models with the desired properties. 

To meet this challenge we construct a general integrable system based on the ancestor Lax operator 
of kunprl and generate in a unified way a series of integrable multi-atom MR models with explicit 
inter-atomic interactions and nondegenerate AF. This includes such new generalizations for JC, BS, 
TI, etc. models and discovers important integrable q-deformations like qBS, qJC, qTI etc. It is worth 
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noting that, our integrable TI model exhibits fuh exponential nonlinearity without any approximation 
and multi-atom qBS and qJC models involves quantum group spin operators. Moreover, since our 
construction is based on a general Yang-Baxter (YB) algebra, together with the generation of various 
models at its different realizations, we can solve them exactly in a unified way through algebraic 
Bethe ansatz (BA). Our strategy of construction is to start with a Lax operator by taking it as a 
combination T(A) = L^{^)Y\f'' Lj{X), with L^{X) linked with the ancestor model of kuriprl and Na- 
number of Lj{X) related to the spin model jcb§the. By construction it must satisfy the YB equation 
R{X — /x)T(A) (8> T{ij,) = (I (X" T(/i))(T(A) (8) I)R{X — /x), with mutually commuting set of conserved 
operators obtained from the expansion r(A) = trT(X) = J2a C'aA" al^a. For standard MR models, as 
we will see below, the Lax operators are rational type linked with the simplest quantum i?-matrix of 
XXX spin chain alja, while for g-deformed models they are trigonometric type related to the i?-matrix 
of xxz chain xxz. We concentrate first on standard MR models and recall that in the rational case 
the 2x2 ancestor Lax operator may be given as 

with operators s satisfying a quadratic algebra 

= 2m+s3 + m~, [s^ = ±s=^, [m^,-]=0. (2) 

The central elements are expressed through arbitrary parameters appearing in as m"*" = 
C1C2, rn~ = c\c2 + C1C2 and as it is easy to see, their different choice reduces to different algebras: 

i) su{u), at = l,m~ = 0, ii) su{l, 1), at = —1, m~ = 0, 
iii) bosonic, at = 0,m~ = —1, iv) canonical, at = = (3) 

and the corresponding limits yield from (Q) the respective Lax operators. In case i), reduces simply 
to the spin Lax operator 

Liw = [ T^' I, (4) 



Sj, A — Sj-\-Cj 



Our Lax operator constructed as above would generate the set of all commuting conserved operators, 
with higher ones containing increasingly higher many-body interactions. The simplest among them 
is Cat^ = s'^ + X^j^" Sj, while the next in the set aCjy^^i, may be defined as the Hamiltonian of our 
unified MR system: 

Hmr = Hd + Hss + Hss, 

Na 

3 

Hss = aY,(s+Sj + S-SI + [c\ + cl)s^S^), 

i 

Hss = aY,({cl + cl)StS^^+c\S-Sj + c^2StSj) (5) 

Here Hss describes matter-radiation, while Hss, matter-matter interactions. Sj, j = 1, 2, . . . , Na stand 
for an array of Na atoms, each with 2s + 1 levels and satisfy the su{2) algebra, s on the other hand 
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signifies a radiation or a vibration mode and satisfies more general algebra ((2j. In @ the radiation 
frequency ujf and the atomic frequencies ujaj,j = 1,2, ...jA'a are defined through inhomogeneous 
parameters of the Lax operator as 

ujf = ^ Wj, Wj = a{ci — C2)cj, uJaj = u!f — Wj + a{c\ + C2) (6) 
j 

Remarkably, the general model © reduces to a new series of integrable multi-atom BS, JC and 
TI models in a unified way at the limits ii), iii) and iv) of For example, case ii) with the choice 

4 = -4 = l,cl = cl^c, (7) 

yields from ^ the model 

Hbs = cofs^ + J2 i^ajS] + a{s+Sj + S-S+)) + a 5](5r5+ - S+Sr), (8) 

j i<j 

which with a bosonic realization of sii(l,l): s+ = \/iV5^,,s^ = h^/N,s^ = N + and the spin- 
s operator S = '^k-, would represent a new integrable multi-atom BS model with inter-atomic 

interactions and different atomic frequencies. Note that at Na = 1, when matter-matter interactions 
vanish and all AF coincide, (jSJ recovers the known model bsbgthe. However we solve below exactly 
through BA the more general case with nonvanishing interatomic couplings and all different AF: LVaj, 
as defined in ©. 

Similarly, a new integrable multi-atom JC model with matter-matter coupling is obtained from the 
same © under reduction iii), consistent with = a, C2 = 0, c} = c, C2 = —a~^ and bosonic realization 
s~ = h, = b\ = b'^b. We do not present here explicit form of this easily derivable Hamiltonian, 
which yields the known model jcb§the only at Na = 1, when interatomic couplings vanish and all AF 
become degenerate. 

We can generate an integrable TI model with interatomic interactions, again from the same MR 
model © at reduction iv), by fixing the parameter values as = — l,c} = c, C2 = c\ = and 
considering consistent realization through canonical variables as = e^*^, = p -\- x. We present 
here only its Na = 1 form by a suitable combination with the other conserved quantity Ci: 

Hti = K - w/)5" + S'^ + a(e-"5+ + e^S") + H^p, (9) 

with Hxp = ^{p"^ + x^) -|- xp, S = ^ J2k^k^ which is a new integrable multi-atom TI model with full 
exponential nonlinearity without approximation. 

For constructing integrable q-deformed MR models the strategy would be the same; only one has 
to start now from the trigonometric type ancestor Lax operator involving q-deformed operators and 
associated with xxz i?-matrix, the explicit form of which is given in kunprl. For simplicity we present 
here only Na = 1 case with the Hamiltonian 

HqMR = Hd + (s^Sg + SgS^) sin a, 

Hd = -icQCOs{aX) +cs\n.{aX), X = {s^- Sg+oj), (10) 

which represent a new class of MR models with belonging to the quantum group Ug{su{2)) and 
Sg to a more general quantum algebra kunprl. It is important to note, that Sg can yield a variety of 
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q-deformed operators, inducing (|T0|) to generate a number of physically relevant q-deformed integrable 
MR models. 

For example, an integrable q-deformed BS model may be constructed from (fT(l|) at cq = 0, by 
realizing Sg through q-oscillator: = ^J[N]qb]j, s~ = bq^[N]q, Sg = + ^, and quantum spin 

operator Sg by using its co-product xxz : = J2j Q '^''<3'^^<j^(^^>i'^^ , = J2j ^j- Note that at 
s = 1, we get an integrable version of an earlier model qbs. 

Similarly the same general model H10|) with choice cq = i,c = 1 and realization = b\^,Sg = 
bg,s^ = N yield a new integrable q-deformation of the JC model, while under reduction cq = i,c = 
and the same realization through canonical operators as for the TI model, it generates an integrable 
q-deformation of the TI model. By taking higher Na values multi-atom integrable variants of all the 
above q-deformed matter-radiation models can be constructed. 

We emphasize that all MR models presented here, similar to their unified construction, allow 
their exact BA solutions also in a unified and almost model-independent way. In BA formalism 
the diagonal entries t(A) = Tii(A) -|- T22(A) produce all conserved operators, while the off-diagonal 
elements T2i(A) = B{X) and Ti2(A) = C(A) act like creation and annihilation operators of pseu- 
doparticles with the M-particle state defined as |M >b= B{Xi) ■ ■ ■ B(Xm)\0 > and the pseudovacuum 
|0 > through C(A)|0 >= 0. The basic idea of algebraic BA alja is to find the eigenvalue solution: 
t{X)\M >b= ^{^,{'^a})\M >B, for which diagonal elements Tii{X),i = 1,2 are pushed through the 
string of S(Aa)'s toward |0 >, using the commutation relations obtainable from the YB equation. 
Considering further the actions Tii(A)|0 >= a(A)|0 >,T22(A)|0 >= /3(A)|0 >, one arrives finally at 
the eigenvalue expression 

M M 

A(A, {A4) = a(A) n /(A - A^) + /3(A) [] /(Aa - A), (11) 

a=l a=l 

where /(A) is defined through the elements of the i?-matrix as for the rational and as 
for the trigonometric case. Expanding A(A, {A^}) in powers of A we obtain the eigenvalues for all 
conserved operators including the Hamiltonian, where the rapidity parameters {A^} involved can be 
determined from the Bethe equations 

Q'(Aa) Yj fjXb - Xg) 

which follow in turn from the requirement of |M >b to be an eigenvector. Returning to our models we 
find that, the major parts in key algebraic BA relations Ullj) and (|12() . described by i?-matrix elements 
/(A), depend actually on the class to which the models belong, rather than on an individual model. 
Therefore, for all standard MR systems including BS, JC and TI models, /(A) is given by its same 
rational form, while for all q-deformed models like qBS, qJC, qTI etc. by its trigonometric expression. 
The only mo del- dependent parts in these equations, expressed through a(A) and /3(A) are determined 
from our general Lax operator construction, which for the rational class using and (@J) is obtained 
as 

a(A) = {c'iiX + r) + c\)Y[{X-s + c,), 

j 

Na 

/3(A) = (c?(A-r)-c^)n(A + s + c,), (13) 
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where r =< 0|s^|0 > depends on particular realization of ((2)) and s = — < 0|S'^|0 > denotes the atomic 
spin. yields easily the needed forms for BS, JC and TI models, at the corresponding choices of 
the parameters like 0, as we have noted above. Similarly, for q-deformed models quantum extension 
of together with the trigonometric form for /(A) have to be considered. For the solution of TI 
and qlT models however one has to adopt a bit different approach close to that of the Toda chain 
falja, since pseudovacuum is difficult to determine for such models. 

For deriving physical consequences from our constructions, we consider integrable two- level multi- 
atom BS model with inter-atomic couplings, by taking spin-i operator Sj = ^aj in ©. Using the 
full strength of the BA method exact solutions for this multi-atom model with arbitrary Na and 
different w^'s can be derived from the same BA relations (|13() for all excitations, by just tuning 

the parameters involved to their required reduction 0. We demonstrate some novel features in Rabi 
oscillation and VRS by using the BA solutions in the Na = 2 atom case of © • The first excited energy 
spectrum Ei = + 2Ai linked with the cubic Bethe equation (|12|) gives three distinct real roots along 
the resonance line ujf = cja(= ^ai = ^a2), resulting to a triplet structure in the VRS with splittings 
in excitation spectrum: Ei = 2.02,3.02,4.02 for uj = 3.02, a = 1. Consequently, the Rabi-oscillation 
becomes involved (see Fig. la) with three transition frequencies. For small detuning 6 = ujf — uJa the 
roots remain real, while at 6 = =f0.30 two of them coincides, collapsing the Rabi-splitting to the usual 
doublet. The excitation spectrum correspondingly reduces to Ei = 2.05,3.81 (degenerate) at sub- 
detuning and to El = 2.22 (degenerate), 3.99 at super-detuning points, reducing the Rabi-oscillations 
to single frequency mode (Fig. lb). Beyond these detuning points two of the roots become complex 
conjugates, leading to an irregular Rabi-oscillation (Fig.lc). 
0.9 
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FIG. 1. Time dependence of transition probability, showing Rabi oscillation at different detuning points: 
at a) resonance, b) degenerate detuning c) further detuning with complex conjugate roots. 

Higher excitations for this model as well as its A^a-atom extension (jSJ can be solved exactly fol- 
lowing the standard BA formalism presented above. We explore now some subtle points and apparent 
controversies regarding the BA solution of BS models, which also have relevance for other integrable 
models. A common belief, though proved only for specific models ndBA, is that the degeneracy con- 
dition (i.e Xa = Aft) for the Bethe states: ~ solved apparently for any integrable 
model. We however find that for our multi-atom BS model at the resonance point 6 = —uJa = 0, the 
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degeneracy condition, which is equivalent to a(Ai) = (3{Xi) = is indeed fulfilled, yielding a nontrivial 
solution Ai = ^{—ujf it 1), which recovers as well the known spectrum for the standard BS model: 
Em = 2M{u!f — Ai) = MEi, at resonance bs. 

Another apparent controversy regarding multi-atom BS models, which is also generic for many 
other models but not emphasized properly in the literature, arises due to the fact that, the dimension 
of the underlying Hilbert space for these models with Na number of two-level atoms has a upper bound 
2^°, for fixed pseudoparticle number. For the standard BS model it is just 2. Therefore the complexity 
of the problem can not increase further for higher excitations with M > Na- For example, in the BS 
model by diagonalizing the Hamiltonian directly one can easily get the exact energy spectrum for 
arbitrary excitation b^s . However, when we try to solve the same problem through Bethe ansatz 
the solution must become increasingly complicated for higher excitations \M >b , since one has to 
find all M Bethe roots {Aa},a = 1,2,. ..,M as solutions to general Bethe equations (|12|) . which is 
impossible analytically! We resolve this problem in an intriguing way by observing that Bethe state 
\M >B and the energy eigenvalues depend in fact not on M number of roots Aq individually, but 
only on some symmetric combinations of them and moreover, the number of these relevant variables 
does not exceed the dimension of the Hilbert space. For the BS model for example, we find them to 
be only two: Xm,Ym and for deriving them explicitly we introduce an equivalent set of symmetric 
Bethe roots through symmetric combinations of the original M roots fut]ire: si = J2a ^a, S2 = J2ab ~ 
AaAf,, . . . ,sm = ria ^a- Combining suitably BAE 1)12(1 and (|11() for the BS model, expressed through 
symmetric roots we arrive at the equations 

EmXm = A+Xm + MYm, EmYm = ^-Ym + MXm, (14) 

where A± = Mojf it ^(5. On the other hand, expressing the Bethe states through two basic states 
of the model we find \M >b= Xm\M,— > +Ym\M — 1,-|- >, i.e. dependant again on the above 
two relevant variables only. It is easy to see that the action of the BS Hamiltonian on this eigenstate 
also reproduces the same relation ()14() . we have derived from the BAE. Fixing cof = 1 for simplicity, 
we find these two variables explicitly through symmetric Bethe roots in the form Xm = Msm and 
Ym = —MsM-'r\5{sM-i + ■ . . + S1 + 1). From the first of the relations (|T1)) we get the energy spectrum 
as Em = + Mkm, where km = and using both these relations derive the simple equation 
Mkm + 5km — M = 0. This quadratic equation is solved easily to yield Em = Mujf it ((5^ + M^)2, 
recovering the known spectrum of the BS model in the general (5 7^ case. Thus through BA we get 
the explicit result for all higher excitations, also analytically, resolving the raised controversy. Similar 
arguments must hold in the corresponding problem for other models. 

Thus we have proposed through general Yang-Baxter algebra a series of new integrable multi-atom 
matter-radiation models including q-deformed models and solved them exactly through Bethe ansatz 
in a unified way. Integrable trapped ion (TI) and q-deformed TI models introduced here are new, 
while q-deformed Jaynes-Cummings (JC) and Buck-Sukumar (BS) models are multi-atom as well 
as integrable extensions of earlier models qbs^qjc. The proposed JC and BS models are nontrivial 
generalizations of well known models jcbethejasbethe, with the inclusion of inter-atomic interactions 
and nondegenerate atomic frequencies. We find that, contrary to the popular belief, the degenerate 
Bethe states do exist in the multi-atom BS models at the resonance point. Multi-radiation modes can 
be included easily in such models preserving integrability. Identifying the models in real systems and 
experimental verification of the related results presented here, especially in many-atom microlasers 
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nQ3, would be an important problem. 
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